


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1993-03 


Thermal stresses of a trimaterial medium ina 
nonuniform temperature field 


Neibert, Michael Joseph 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/27308 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


/ (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

; | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 
































ot 


‘ ' 
. 
ene 
' 
. 
‘ 
' 
. ' 
' 
* @e 
’ 
1 
. ete 
' : . 
1 
s . 
et 
. 
' 
' 
- € 
a 
. . rae 
1 bor ’ 
oo 
. ' 
' 
oe 
.e 
tay 
Ae 
' 1] 
' 
' 
. 
' 
‘ 
a 
’ ? 
. 
' 
° 
' 
i . 
. 
‘ . 
' 
' 
' 
hd 
. 
' 
) 
oe 
. 
' 
, . 
' ' 
F ‘ 
° 
' 
' ' 
“4 
1 1 
Jar 
Ld - 
' 
et 
° 
° 
ae 


“—s 





7 








a) 
\ en 
Sete tiggh O] 
' atee 
ia Pee ‘ uly fired 





. 




















aie) Pe Rak: 
oY rey 1 





" seit i 
hebetdh 45 aAsaibetane Ose B 
4‘ BOE Fe est h Aho ok ay #4 
of owas Bae Ae ye ott grt ‘. Peden 740 
iateun PG ie iaod dads betaned ote 
puatere todos: ee? Coors 
Geel oo ¥ 

gy tet Oph BiuSe ateee! 5.0, My eri 
Beat a 4 ? batt he anit ie 
ware ay ea ba ibat we 

Pah dicate astartles« aa ie 

Are ore be 






RA 































































tesa aie ae & 

y) ea siatisired 4 wed § 
AL nat Jy ae ivy a 7s vitae i « £y.8 Tm wantet ttn 
ee ee porta LIN : 
va! oF} ba Meet ltae ih 

u 


4 











tae 

194 Saiit dat 
i 

wai, ae ry oh 





va 
23 
eam en 





o* asa? 
ath@. 






pan atea ) fd 
a 





SAN ot tne dean 
“tte atet 
“Diet as Be ies * 
safe ef sue te 
ty RL gM4 fel 
seali. aw whe ae 


tote. Cate 


ee oe ee Decent 
e808 obbet ine 

















er Ord te 
ohehed, ake 




















ai « 
*oatek ye ot 
itis: Sa 
wAmshio! 
v tyr. Satesstanctors " 


hase oo: Seder 















Fosses a bedtee I 
saes.e, prerihe = 
taredarle No ‘cag: 















. 


Wweae wt as 
. 





r “I 
ry 





BF 

















Te. 
wr 


. 
= 
~aa Se 


= a 




















































































2 ate yey 
e P Cia ms AG F] 
' . . 
ats . ¢ . ’ . " oe 
1 ‘ Pee ieee ies tn : 
aca ” . . 
a oe é ' t ’ 
' LC . . 
ne Ce | 8 ie i ’ 
te. spehaety Grea ue "hee 
' eae : ae fs A . tity meen gtt g 
a ‘ os : eta Se8 
. . . 
: J ’ fa ef 1 . ‘ oe 
pute Bae « 1 Pan oan 
’ . ' ee ' 
‘ 
1 . of 8 . 1? a om . 
. ‘8 . ‘ 
oe ' . : . 
a ' ; migaoe : La | : 
eo Ce ee er ‘ ’ a 
Tot A Ay Pes : Tear 
a a! i ‘ . ‘ 
' ' ’ ‘ 
er : ii ley 
« eo. - . . etie st 
' ' s . 
: ‘ ti . it ” - . 
’ r ‘ 1 Fi 1 . a « »* 
. i ie Be a OF : By fasts 
" ' . LORY Te 
aan ® . . ‘ 
aa oe. oe 5 
' * EIT i 
t : de (every ae 
' : meee : ; ise ‘< “2% it Ras 
' 1‘ ee ‘ ' ‘wa 
“e @ i aseie ig ye Vee este 
’ ' oan ‘ "Su 
of ’ 8 owe 
*s > : pe 7 at 
2 i ar . e oh Sy Slat egeee tl 
: ‘ ne =< ; tr nee) te 
S. hs “5 ‘ F Dae « cae Ned : 
. ' ® . 4 i‘ 
: . Py ‘ ' rie ' ( “ ‘ so tglpiucnn shaban bea 
Pe en er Roa T aspen ct 
1 ae « . - a 
- at ‘ Sg Oe ee Ane aes ae ' ; : s * Set reagan eee 
' iS . wee oe . « 
. “. ' 6) ce ve : 
¢2 
nae oy ' chree see, Wightace 
aay ere ae ae Pet te PNR Rant 
‘ 
a | . Caan fee: 
u ak tia” sole 
sg Phan ’ Pla Tin eo ty "33 tare 
Y . ' hs ' be aes 
oh UJ o be ope ' eh: 
oat ls ite . A Sav G9 ay A 
S See 8 ba 9 x 7 Witery Saivtores 
‘ - ' ve $d sine Ae Fae Red ie ta 
‘ a al 
ac free : : 7 a) 7 sae 
= 1 ° ° ae ' 
oe ‘ 1 oe aac ‘ : be 
sey Tn es oo arte 
. * . . 
. ' $ a 
7 ' oe i. a . A ate = 
' pe 1 ae ' oe 
¢ ' 
rer) " : Miler ae 
+ to. 
ome . . ? a 1! ‘ ry cal gis eghins i vs " ooeses 
ro oo Si OE! va of tae Bus * 
; ’ . at - ra 
: (9 ‘8 a a | ~ « A a stn te 
2 nee ‘ BS 4 ates s 3s nip 
‘ 4 Pia ( " gae ie eres relia sh weaeee = 
i : . ste ate ie 
Seth * : Saar gtiiniats t, tiga eete ty p97 KY Fey sei “ * 
‘ 2 . a i we pert beh fe Bt efron By : B4 Bp ut 
eae Pag Oe eb ip) wre: cig’ S45 twSy pee es i mae see: 
. ne 2 tat! 
' 7 - . ‘| 
4 ‘ ’ ' of 
. ‘ 7 ‘ 
‘ . aS ‘ zn a ‘ 
. : - 1 . ’ Penees ie 
a 7 4 ‘6 ert 4 | 
. Gi oy E ty a 
, ' at ' ‘ e) moat 
' 2th 8 A Z oan 
: , see : al pghiaa ty py apse tees a 
' 5 ; F : i¢ wpe 409. sete Wan #9? ritP- 
. ot aoe id 4 Tas oy 5 
2 r e . rs 4 Pe ye s ue 





DUDL=VY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY CA 93943-5101 














Approved for public release: Distribution is unlimited. 


Thermal Stresses of a Trimaterial Medium in a Nonuniform 
Temperature Field 


by 
Michael J seis ar 
Lieutenant, United States Navy 
B.S., Colorado State University, 1981 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN MECHANICAL ENGINEERING 
from the 


NAVAL POSTGRADUATE SCHOOL 
March 1993 


Form roved 
REPORT DOCUMENTATION PAGE 


Public reporting burden for this collection of information is estimated tO average 1 hour per response, including the time for reviewing instructions, searching existing data sources 
gathering and maintaining the data needed, and completing and reviewing the collection of information Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and 8udget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 


1. AGENCY USE ONLY (Leave blank) | 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 
March 1993 Master's Thesis March 1993 


4. TITLE AND SUBTITLE 5. FUNDING NUMBERS 





Thermal Stresses of a Trimaterial Medium in a Nonuniform 
Temperature Field 


6. AUTHOR(S) 
Michael Joseph Neibert 





7. PERFORMING ORGANIZATION NAME(S) AND ADORESS(ES) 8. PERFORMING ORGANIZATION 


REPORY NUMBER 






Naval Postgraduate School 
Monterey, CA 93943-5000 





10. SPONSORING / MONITORING 
AGENCY REPORT NUMBER 


9. SPONSORING /MONITORING AGENCY NAME(S) AND ADDRESS(ES) 





Naval Postgraduate School 
Monterey, CA 93943-5000 





11. SUPPLEMENTARY NOTES 
The views expressed are those of the author and do not reflect the official policy or position 
of the Department of Defense or the U.S. Government 





12a. DISTRIBUTION / AVAILABILITY STATEMENT | 1 12b. DISTRIBUTION CODE 


Approved for public release: distribution is unlimited 








13. ABSTRACT (Maximum 200 words) 


The objective of this investigation was to conduct a parametric study of the effect of a 
nonuniform temperature field on the system behavior of a trilayered medium. In particular, 
the resulting shear and normal stresses along the media interfaces are analyzed. A finite 
element model utilizing a recently developed element which provides for both axial and 
lateral displacement continuity is employed. First, the effect of the material properties, that 
is, Young's Modulus and coefficient of thermal expansion, is examined. Then, the effect of 
the geometric properties, that is, length and thickness dimensions of the midlayer, is 
analyzed. Finally, a study of the effect of a nonuniform temperature field on the trimaterial 
medium is conducted. 


. SUBJECT TERMS 15. NUMBER OF PAGES 
Thermoelastic Stresses, Multilayered Media, 108 


Nonuniform Temperature Fields 16. PRICE CODE 





19. SECURITY CLASSIFICATION | 20. LIMITATION OF ABSTRACT 


OF ABSTRACT 
unclassified 


NSN 7540-01-280-5500 ° Standard Form 298 (Rev 2-89) 
1 Prescribed by ANSI Std 239-18 


SECURITY CLASSIFICATION 
OF REPORT 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 






unclassified unclassified SAR 





ABSTRACT 


The objective of this investigation was to conduct a parametric study of 
the effect of a nonuniform temperature field on the system behavior of a 
trilayered medium. In particular, the resulting shear and normal stresses 
along the media interfaces are analyzed. A finite element model utilizing a 
recently developed element which provides for both axial and lateral 
displacement continuity is employed. First, the effect of the material 
properties, that is, Young's Modulus and coefficient of thermal expansion, is 
examined. Then, the effect of the geometric properties, that is, length and 
thickness dimensions of the midlayer, is analyzed. Finally, a study of the 
effect of a nonuniform temperature field on the trimaterial medium is 


conducted. 


II. 


IIT. 


TABLE OF CONTENTS 


INTRODUCTION .&. ..c . = ee 1 
FEM MODEL FORMULATION ................... 733 +) 
A. PROGRAM DEVELOPMENT ................... 7 oo 5) 
B. PROGRAM INPUT .............9. 999399 12 
C. PROGRAM OUTPUT ...s...... . . SRR 13 
EFFECT OF MATERIAL PROPERTIES ..................... 14 
A. CASES WAND 2 .. ww... es. 16 
B. CASES 3, 4, AND 5 0.60... ...........-.--+: 1 or 20 
C. CASES 3, 6,7, AND 8 ......................., —o—oo 24 
D. CASES 7,9, 10, AND 11 .....%............°07 7 ooo 29 
KE. SUMMARY ............ ee eee eee ee 29 
EFFECT OF MIDLAYER THICKNESS ...................... 36 
A. NORMAL STRESS DISTRIBUTION ..................... 36 
B. SHEAR STRESS DISTRIBUTION ................ 733m 39 
C. SUMMARY .....5.....004 52555-5159 45 
EFFECT OF MIDLAYER LENGTH ................/ 49 
A. NORMAL STRESS DISTRIBUTION ............... 7 49 
B. SHEAR STRESS DISTRIBUTION ................. 33 52 
C. SUMEMARY . woo. 6 eb ee nt ee sg oe 56 


1V 


DUDLEY KNOX LIBRARY 


NAVAL POSTGRADUATE 
SCHOG 
MONTEREY CA 93943.5; ts 


VI. EFFECT OF NONUNIFORM TEMPERATURE DISTRIBUTION .. 57 


Pon veiy STRMoomISURIBUTION ...................-- 08 
Beou hake iRmeeeISTRIBUTION ..................200.. 60 
a 64 
WIT SUMMARY AND CONCLUSIONS ©... 2... ee 66 
A. EFFECT OF MATERIAL PROPERTIES ................... 66 
B. EFFECT OF MIDLAYER THICKNESS ................... 67 
PeePPECT OF MIDLAYER LENGTH ...................... 67 
D. EFFECT OF NONUNIFORM TEMPERATURE DISTRIBUTION 68 
E. RECOMMENDATIONS FOR FURTHER RESEARCH ........ 68 
Pure OIX AlN PROGRAM 2... ee 69 
PuerNDIX By PROGRAM EXEC 2... 1... ee ee 76 
PemeNOIX CC; INPUT PIE o.oo ee ec eee 77 
ee IBID GIB) (Ce ie G1 ial) b,c 81 
Peete REMORMINGHS ...... 0M. cee LOL... 98 
Mee WMS TRIBUTION LIST «ww... eee g9 


LIST OF FIGURES 


Figure 1.1 Electronic Package ................... se 2 
Figure 1.2 Tri-Material Configuration .......................... 4 
Figure 2.1 A Typical Element With Six Degrees of Freedom ......... 5 
Figure 2.2 Finite Element Analysis Mesh for Tri-Material 

Configuration .............+..+..+.sssss, ) 13 
Figure 3.1 Tri-Material Configuration ..................,..— 14 
Figure 3.2 Normal Stress at Interfaces for Cases land2 .......... lg 
Figure 3.3 Shear Stress at Interfaces for Cases 1 and2............ 18 
Figure 3.4 Shear and Normal Stresses on the Midlayer for 

Cases land2...........000...0055... 9 
Figure 3.5 Normal Stress at Interfaces for Cases 3, 4,and5 ........ 21 
Figure 3.6 Shear Stress at Interfaces for Cases 3, 4,and5 ......... 22 
Figure 3.7 Shear and Normal Stresses on the Midlayer for 

Cases 3, 4,and5 .................2++++0+. 23 
Figure 3.8 Normal Stress at Interfaces for Cases 3,6, 7,and8 ...... 25 
Figure 3.9 Shear Stress at Interfaces for Cases 3,6, 7,and8 ....... 26 
Figure 3.10 Shear and Normal Stresses on the Midlayer for 

Cases 3 and 8 1... .... cc ete tee + er ae 
Figure 3.11 Shear and Normal Stresses on the Midlayer for 

CasesGand7........552...20022+++2+++ ) 28 
Figure 3.12 Normal Stress at interfaces for Cases 7,9, 10, and 11 .... 30 
Figure 3.13 Shear Stress at interfaces for Cases 7,9, 10, andll...... 31 


Figure 3.14 


Figure 3.15 


Figure 4.1 


Figure 4.2 


Figure 4.3 


Figure 4.4 


Figure 4.5 


Figure 4.6 


Figure 4.7 


Figure 4.8 


Figure 4.9 


Figure 5.1 


Shear and Normal Stresses on the Midlayer for 
ee eee. ©... Sy. a oes. 32 


Shear and Normal Stresses on the Midlayer for 
acc cmi@iamee!) lwweeie . eee. erm ee ee oo 


Tri-Material Configuration for Cases 5, 12, 18, 14, and 15 . 37 


Normal Stresses at the Interfaces for 
Oe sg 1 Ie ee 38 


2-D and 3-D Plots of Normal Stress/Maximum 
Normal Stress at the Upper Interface for 
Cecee 20M (bee sit ha 40 


2-D and 3-D Plots of Normal Stress/Maximum 
Normal Stress at the Lower Interface for 
Caco mee eel Among oe fe... ws. we ew 41 


Shear and Normal Stresses on the Midlayer for 
Caseses, (Zreameeihmeren., . WP. . Pe ee 42 


Shear and Normal Stresses on the Midlayer for 
Cae aR ayich 1 Ga 43 


Shear Stresses at the Interfaces for 
Waceotoe cero laeaad 1D... kk te tt ht tt ee 44 


2-D and 3-D Plots of Shear Stress/Maximum Shear Stress at 
the Upper Interface for Cases 5, 12, 18, 14, and 15....... 46 


2-D and 3-D Plots of Shear Stress/Maximum Shear Stress at 
the Lower Interface for Cases 5, 12, 18, 14,and15....... 47 


Tri-Material Configurations for 
Cegee 1% We USMS) Eta] 710 50 


Figure 5.2 Normal Stresses at Interfaces as a Function 


Figure 5.3 


of Non-dimensional Location for 
Saeccel 2 IG illo. WoweamenZO ww we ee we eee 51 


Shear and Normal Stresses on the Midlayer for 
oe siee. ei, Birael )70 g 53 


Figure 5.4 


Figure 5.5 


Figure 6.1 
Figure 6.2 


Figure 6.3 


Figure 6.4 


Figure 6.5 


Figure 6.6 


Shear and Normal Stresses on the Midlayer for 


Cases 18, 19, and 20 .......... ee.) 54 
Shear Stresses at Interfaces as a function of Non-dimensional 

Location for Cases 12, 16, 17, ‘18, 19, and 20............ 55 
Tri-Material Configuration for Cases 21, 22, 23, and 24 .... 57 
Material A Linear Temperature Distribution for Case 24 ... 58 


Normal Stresses at the Interfaces for 
Cases 21, 22, 238%and 24 .°9...... 3) ae ee 59 


Shear and Normal Stresses on the Midlayer for 
Cases 2lmande22 ........... 088 .... 08... | 61 


Shear and Normal Stresses on the Midlayer for 
Cases 23 andi24 ..................... 5 62 


Shear Stresses at the Interfaces for 
Cases 21, 22, 28, and 24 .............-...-+.))3——— 63 


I. INTRODUCTION 


The history of thermoelastic stresses on layered materials begins with 
Timoshenko. In a 1925 paper [Ref. 1], he first discussed the case of a directly 
bonded bi-material configuration. He studied bi-material thermostats modeled 
in the form of high aspect ratio beams subjected to uniform heating or cooling. 
Bearing stresses in closed form expressions were presented and it was noted 
that there was not an elementary way of determining the shear stress 
distribution along the bearing surface. The curvature of the beams at a given 
temperature was shown to be related to the beam geometry, the elastic 
constants of the materials, and their coefficients of thermal expansion. 

Goland and Reissner [Ref. 2] and Christensen [Ref. 3] studied the shear 
and peeling stresses along the interface of single lap adhesive joints subjected 
to axial loads. Goland and Reissner derived closed form analytical solutions 
which showed that the peak peeling and shear stresses occurred at the joint 
edges. Additionally, it was shown that all stresses depended on the joint 
dimensions and the stiffness properties of the plates. Christensen performed 
his analysis using the finite element method. He determined that a large 
stress gradient existed across the thickness of the adhesive as well as along 
the lap, especially near the ends of the lap. Neither study investigated the 


effect of temperature on system behavior. Burgreen [Ref. 4] analytically 


computed interfacial shear stresses in layered beams and examined the 
behavior of a bi-metallic beam subjected to a transverse temperature variation. 
Lau [Ref. 5] derived closed form stress analysis solutions for a chip in a semi- 
infinite substrate and a chip on substrate with finite thickness subjected to 
temperature rises. 

A typical electronic package consists of a semi-conductor (chip) attached 
to a substrate material by an adhesive as depicted in Figure 1.1. Placing more 
power in an ever decreasing space is the trend in electronic component 
packaging. Power densities of 12.5 MW/m° have been reached by today's 
integrated circuit technology [Ref. 6] and future integrated circuits are 
expected to have at least twice this power density leading to higher operating 
temperatures. The resulting high operating temperatures give rise to thermal 
cycling causing thermal fatigue at the different material interfaces and can 
lead to premature failure of the aa To prevent failure, a maximum 
temperature of between 100°C and 110°C at a semiconductor junction has been 


set by the military [Ref. 7]. 
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Figure 1.1 Electronic Package 
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While conducting a literature search, numerous studies concerning the 
topic of thermal response of electronic packages was uncovered. In perhaps one 
of the most extensive studies, Suhir [Ref. 8] reviewed work conducted on bi- 
material and tri-material electronic packages that were either subjected to a 
uniform temperature or an applied axial load. For the bi-material configura- 
tion, it was determined that the maximum peeling and shear stresses along 
the interface occurred at the outer edges of the assembly. Additionally, it was 
shown that the stresses depended upon the axial and interfacial compliances 
of the component layers. 

Royce [Ref. 9], Riches [Ref. 10], and Chung [Ref. 11] all state that the 
primary cause of mechanical stresses that are induced in electronic packages 
are a result of the difference in coefficients of thermal expansion of the 
different materials in the assembly. Royce commented that thermal cycling 
during normal operation could give rise to thermal fatigue at the various 
material interfaces and could lead to premature failure of the structure. Rich 
reported that the types of stresses that arose from thermal mismatch consisted 
of stresses normal to the chip and interfacial stresses. He predicted that 
normal stresses in the chip were at a maximum in the chip center and 
decreased to zero at the chip edge. However, the interfacial stresses were zero 
at the center and rose to a maximum at the edge. Rich concluded that 
compliance and thickness of the chip attachment materials were critical to 


controlling the stresses imposed. Two ways to reduce or avoid thermally 


induced stresses in electronic packages was proposed by Chung. The first was 
stress reduction by coefficient of thermal expansion matching and the second 
was stress reduction by the use of extremely flexible adhesives. 

The purpose of this investigation was to gain further insight and under- 
standing into the effects that the various system parameters have on thermally 
induced behavior of an electronic package as shown on Figure 1.1. 

To achieve this, the electronic package was generalized into a tri-material 
configuration as displayed in Figure 1.2 and a parametric study was conducted 
using finite element analysis. First, the effect of the material properties, 
Young's Modulus and coefficient of thermal expansion, was examined. Next, 
the length of the midlayer was varied. Then, the effect of midlayer thickness 
on system behavior was studied. Finally, nonuniform temperature fields were 


imposed on the tri-material configuration. 
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Figure 1.2 Tri-Material Configuration 


Il. FEM MODEL FORMULATION 


The purpose of this investigation was to determine the stresses resulting 
from a nonuniform temperature field on a tri-material configuration. To 
accomplish this, a finite element method (FEM) program was constructed. 
This program is presented in Appendix A and the program exec is in 
Appendix B. 

A. PROGRAM DEVELOPMENT 

A recently developed element which provides for both axial and lateral 
displacement continuity was utilized in the stress formulation. As shown in 
Figure 2.1, each element has six degrees of freedom; axial displacements at the 


four corner points; and lateral displacements at the two ends. 





Figure 2.1 A Typical Element With Six Degrees Of Freedom 


The lateral displacement field is determined by linear interpolation 


V°=N,V; tN,V> (2.1) 


where v, and v, are the lateral displacements at the ends and where the linear 


shape functions N, and N, are given by 
N,==(-a) N,=# (2.2) 


resulting in 


2 
v(x)=)> Na); (2.3) 
i=l 
The axial displacement field is assumed to be linear in both the axial and 
transverse direction. That 1s, 


u " =H i hg +H iu : (2 ° 4) 


iis =Hy uy +Hius (2.5) 
where superscripts b and t refer to the bottom and top respectfully. The linear 


shape functions H, and H, are given by 
Hp=Hy=—(h-y) Hj =Hy=4 (2.6) 


resulting in 


2 
u(xy)=)> N,@)LH, (yu; +H(0)u;) (2.7) 


i=l 
From strain-displacement relations 
Ou oO 


aa EN, (x) Hu; SDN, (x)H,(y)u; (2.8) 





or 


: ON; b t 
ora = K; 2 u;) 
2, ax l Zz 


and 
_ Ou. OV 
ies A 
oy ox 
Hy AH, u, Hyu;) +— 
oye tod SP nn} 
or 


2 OH, OH 2 ON. 
Yay SNe : Sauls me 


i=l 


The derivatives of the linear shape functions are 











ON, 1 ON, 1 OH, 1 OH, 1 
ae | a~tlwyh wh 
This results in 
H,(y) _ 
: (hy -u)+——(u,-u)) 
and 
N,@) mi va 
yF tu; - u, + {u5- uy}+ 


(2.10) 


(211) 


(Ao) 


(2.14) 


(2915) 


Defining the axial displacement vector as 


{8 )7=Uy uy uy us) (2.16) 
and the stiffness matrix as 
lph 
e Zl, 
[K,] f, I {BIE(BY dydx (2.17) 
where the B vector is 
ON ON ON. ON. 
(Byias(—_§$'H —'9. —42H —2H (2.18) 
ae | atm mw tO Ok 


and the force vector due to temperature as 
_Eh eich 
Fy ff, BiEa A Tdydx (2.19) 


gives the bending matrix equation as 
[K,]{6 J =F, (2.20) 


Equation 2.20 defines the bending behavior. Additionally, the stiffness matrix 


can be shown as 


2 1 @.-1 
1 2 -1 -2 
ae (220) 
6l)/-2 -1 2 1 
1 -2 1 2 


Behavior due to shear is obtained as follows. Define the row vector of 


displacement degrees of freedom as 


8) =U) Uy V) uy Uy v,) (2.22) 


The shear stiffness matrix is given by 


[K3=[ f (BIGB dadydx (2.23) 


where 


OH OH ON OH OH ON 
ee NS ON OC) (2.24) 
oy dy = o& oy dy o& 
which gives the equations for shear behavior 
[K,]{6 }=(0} (2.25) 
Additionally, the stiffness matrix can be shown as 
1 
4h 4h 2 4h 4h 2 
a, 2 ee 
4h 4h 2 4h 4h 2 
tort A 1 al sh 
[K]=G 2) ie Sia) as | (2.26) 
ae el gl 
4h 4h 2 4h 4h 2 
aot at ti ol 
4h 4h 2 4h 4h 2 
mee yet 1 1 ok 
eae 2 2 
The system matrix for bending and shear is 
[A]=[K,] +[K;] (2.27) 


Let 


then 


-d dad -b -d “Wed 


The force vector due to a temperature field is 
sae ne 
F4=0°f [ (Bldyde 
where 


Q°=(Ea AT) 


or 


(F4}=0° dydx 


10 


(Zize) 


(2.29) 


(2.30) 


(2.31) 


(2732) 


Let 


Ac=ae! (2.33) 
2 
then 
-] 
-] 
0 9 
iF 4=/A® (2.34) 
1 
0 
After solving 
[K]{S}=tF} (2.35) 


for {5}, the stresses can be calculated. For the e” element, the bending stress 








is 
g,=Ee. (2.36) 
or 
ay. 
see canes! (2.37) 
l l 
my 5-6 
of=— (uy-u))= - 2 (2.38) 


where 6,” and o,” are the bending stresses acting on the e” element at the 
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bottom and top respectfully. The shear stress is 








me eee (2.39) 
or 
b 
Ro a aL Md (2.40) 
h i h [ 
t b 
ot Gla Be esees “amas (2.41) 
2 h l h l 


where coe and Te are the shear stresses acting on the e“ element at the 
bottom and top respectfully. Finally 
7 § Sass P ae (2.42) 

B. PROGRAM INPUT 

The input file for the finite element program is displayed in Appendix C. 
A mesh grid with 200 elements as shown in Figure 2.2 is generated. The input 
file allows for variation in both the element layer thicknesses and the element 
row widths. Values for Young's Modulus, Poisson's Ratio, and the coefficient 
of thermal expansion can be inputted row by row. Additionally, the 
temperature for each level of the mesh can be inputted. Appendix C is 
annotated, showing the placement of the different inputs. All dimensions are 


in millimeters. 
12 


C. PROGRAM OUTPUT 

The output file from the finite element program is shown in Appendix D. 
The force and displacement vectors for each node are outputted. The top and 
bottom bending stresses and shear stress for each element are also displayed. 


Additionally, the normal stresses along the material A and material B 


interface and the material B and material C interface are outputted. 





| 
pee tHe 


Figure 2.2 Finite Element Analysis Mesh For Tri-material Configuration 
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Il. EFFECT OF MATERIAL PROPERTIES 


This chapter studies the effect of Young's Modulus and the coefficient of 
thermal expansion on the tri-material configuration of Figure 3.1. Several 
computer runs were made assigning different Young's Modulus values to the 
three materials while maintaining a constant coefficient of thermal expansion 
value. No stresses resulted. As a consequence, the remainder of this chapter 
is concerned with variations in the coefficient of thermal expansion (a). Table 
3.1 shows the parameters used for the 11 cases that were evaluated. For all 
cases, Young's Modulus was given a constant value of 100 GPa for the three 


materials and a change in temperature of 100°C was imposed. 


“MATERIAL. 


MATERIAL C 





ALL DIMENSIONS IN MILLIMETERS 
CENTERLINE 
Figure 3.1 Tri-matenial Configuration 
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TABLE 3.1 COEFFICIENT OF THERMAL EXPANSION STUDY 


| 4 | 00 | 50 | 00 








The eleven cases were grouped into four combinations to allow 
comparisons to be made. The combinations were: 


e Cases 1 and 2. a, and G&, were set equal and a, was allowed to range 
from one-half a, &_ to twice Oa, Qo. 


e Cases 3, 4, and 5. ,, Gp, and G were given different values in different 
combinations. 


¢ Cases 3, 6, 7, and 8. O, and a, were held constant but at different values 
while a, was allowed to range at values greater than a, and M. 


¢ Cases 7, 9, 10, and 11. a, and a, were held constant but at different 
values while a, was allowed to range at values less than a, andag. 
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A. CASES 1 AND 2 

Cases 1 and 2 investigated how the tri-material configuration responded 
when the top and bottom layers (materials A and C) held the same value of the 
coefficient of thermal expansion while the middle layer was given values lower 
(case 1) and higher (case 2). 

Figure 3.2 is a plot of the normal stress along the upper interface (between , 
materials A and B) and lower interface (between materials B and C) as a 
function of position. For case 1, both the upper and lower interfaces are in 
tension, showing that the top and bottom layers with their larger amount of 
expansion, wanting to peel away from the middle layer. In case 2, the normal 
stresses at the upper and lower interfaces are all compressive, that is, have 
bearing stresses. The middle layer wants to expand more than the top and 
bottom layers and pushes against them. For both cases, the maximum stresses 
occur at the ends of the configuration. 

Figure 3.3 is a graph of the shear stresses acting on the upper and lower 
interfaces. The shear stresses for case 1 tend to axially extend the midlayer 
while the shear stresses for case 2 tend to axially compress the midlayer. 
Figure 3.4 shows the shear and normal stresses acting on the midlayer for 


cases 1 and 2. 
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B. CASES 3, 4, AND 5 

Cases 3, 4, and 5 investigated the response of the tri-material configura- 
tion when all three layers had different values for the coefficients of thermal 
expansion. Case 3 was arranged with magnitudes as: &,> Q,> Qo, case 4 as: 
O,> Ao> Ap, and case 5 as: p> A,> Ac. 

Figure 3.5 is a plot of the normal stresses as a function of position along 
the upper and lower interfaces. Case 3, with the coefficients of thermal 
expansion in decreasing order from top to bottom, displayed minimal peeling 
stresses along the interfaces. Case 4, with the coefficient of thermal expansion 
lowest in the midlayer, exhibited characteristics similar to case 1. The normal 
stresses were peeling on both interfaces. Case 5, with the coefficient of 
thermal expansion in the middle layer the largest, behaved like case 2. The 
normal stresses along the upper and lower interfaces were bearing throughout. 

The plots for the shear stresses along the upper and lower interfaces are 
shown in Figure 3.6. In case 3 the shearing stresses were the lowest, with 
material B wanting to elongate axially along the upper interface. The 
midlayer wanted to axially elongate for case 4 and to axially compress for case 
5. This was the same behavior that was observed for cases 1 and 2, respectful- 
ly. Figure 3.7 displays the shear and normal stresses acting along the 


midlayer interface. 
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C. CASES 3, 6, 7, AND 8 

This study compared the responses of the tri-material configuration when 
OQ, and &, were held constant at values lower than a, (with > O,), and @, was 
then allowed to increase as follows: a, (case 6) < a, (case 7) < O, (case 3) < a, 
(case 8). 

Figure 3.8 is a graph of the normal stresses along the upper and lower 
interfaces for the midlayer. When the value of the coefficient of thermal 
expansion for the upper layer was kept reasonably close to the coefficient of 
thermal expansion of the midlayer (within a factor of two) as shown by cases 
6 and 7, the normal stresses tended to be minimal and were mainly bearing. 
However, when a, was further increased as in cases 3 and 8, the magnitudes 
of the normal stresses increased and became peeling throughout the upper and 
lower midlayer interfaces. 

The shear stress plot for the upper and lower interfaces is given in Figure 
3.9. Cases 6 and 7 resulted in shearing stresses at both ends of the upper and 
lower interfaces with minimal shear stresses along the majority of the 
interface lengths. In cases 3 and 8, the shearing stresses on both interfaces 
tended to axially stretch the midlayer with case 8 having a larger magnitude 
of stress. Figures 3.10 and 3.11 display the shear and normal stresses acting 


on the midlayer for cases 3, 6, 7, and 8. 
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D. CASES 7, 9, 10, AND 11 

In this grouping, the coefficients of thermal expansion for materials A and 
B were held constant at values greater than material C (with a, > a,), and a, 
was allowed to increase as follows: OG, (case 9) < A, (case 7) < a. (case 10) < a 
(case 11). 

Figure 3.12 is a plot of the upper and lower interface normal stresses for 
the four cases. Similar characteristics were displayed by cases 7 and 9. Along 
the upper interface the normal stress was initially peeling, decreased to zero, 
and then became bearing. At the lower interface the normal stresses were 
initially bearing, went to zero, and then became peeling. For cases 10 and 11 
the normal stresses were peeling on both the upper and lower interfaces. 

A plot of the shear stresses acting on the midlayer is shown in Figure 
3.18. As can be seen, along the upper interface all four cases displayed similar 
shear stress characteristics with the surface being pulled in axial tension. 
However, at the lower interface, cases 7 and 9 showed axial compression while 
cases 10 and 11 displayed axial tension. The shear and normal stresses acting 
on the midlayer for these four cases are shown in Figure 3.14 and 3.15. 

E. SUMMARY 

A review of the peeling and shear stresses acting along the upper and 

lower interfaces reveals force equilibrium in both the lateral and axial 


directions. However, moment equilibrium is not satisfied exactly. This 
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deviation is most probably a result of the type of finite element that was 
chosen. 

As expected, when materials A and C were given the same properties, 
symmetrical stress conditions resulted. This was clearly shown in cases 1 
and 2. When the coefficient of thermal expansion in the midlayer was smaller, 
the upper and lower layers wanted to separate from the middle layer. When 
the coefficient of thermal expansion for the midlayer was greater, the top and 
bottom layers tended to compress it. 

Maximum peeling stresses occurred when the coefficient of thermal 
expansion for the top or bottom layer was much greater than the middle layer 
as shown in case 4. Similarly, the maximum bearing stresses arose when the 
midlayer coefficient of thermal expansion was greater than the upper and 
lower layers as evidenced by case 5. 

The relative difference in the coefficient of thermal expansion between the 
layers played a similar role for the shear stresses. The maximum tendency to 
axially elongate the midlayer occurred when the coefficient of thermal 
expansion for the middle layer was much less than the upper or lower layers. 
Additionally, the maximum tendency to axially compress the midlayer 
happened when the coefficient of thermal expansion for the midlayer was much 


greater than the top or bottom layers. 
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Minimal normal stresses resulted when the values of the coefficients of 
thermal expansion for the upper, middle, and lower approached each other. 
The conditions for minimum shearing stresses were the same. 

The study of cases 7, 9, 10, and 11 produced an interesting result. As can 
be seen, cases 7 and 9 exhibited very similar stress states even though the 
coefficients of thermal expansion for the bottom layer differed by a factor of 
ten. More computer runs were conducted, further reducing the coefficient of 
expansion by several factors of ten. The changes in the state of stress was 
minimal. 

Finally, when the coefficient of thermal expansion for the tri-material 
configuration was kept constant and uniform, changing of the Young's Modulus 


values did not produce any resulting thermally induced stresses. 
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IV. EFFECT OF MIDLAYER THICKNESS 


This chapter investigates the effect of changing the thickness of the 
midlayer (material B) on the tri-material configuration. Five different 
midlayer thickness cases were studied and compared. As displayed in Figure 
4.1, the following thicknesses for material B were used: 1.00mm (case 5), 
0.50mm (case 12), 0.25mm (case 13), 0.10mm (case 14), and 0.05mm (case 15). 
For all cases, the dimensions of materials A and C were kept constant (10mm 
length and 1.0mm thickness). Additionally, the length of the midlayer was 
held at 5mm and the tri-material configuration was given a temperature 
increase of 100°C. The material properties for the five cases were given the 
following constant values: E, = Eg = Eco = 100 GPa, o,= 100x10°/°C, ap= 
300x10°/°C, and a = 10x10°/°C. 

A. NORMAL STRESS DISTRIBUTION 

Plots of the normal stresses along the upper and lower interfaces of the 
midlayer for the five cases are displayed in Figure 4.2. Both the upper and 
lower normal stress distributions exhibited similar characteristics. The 
stresses were always bearing along the interfaces for all five cases with the 
maximum stresses occurring at the left and right ends of the midlayer. For 
every case, the minimum bearing stress was developed midspan of the 


interface. As the thickness of material B was decreased, the bearing stresses 


36 








Figure 4.1 Tri-Material Configuration For Cases 5, 12, 18, 14, And 15 
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similarly decreased. However, the shape of the normal stress distribution 
curves remained comparable. 

Comparisons were made between the maximum of the normal stresses of 
all five cases on the upper and lower interfaces and the remainder of the 
normal stresses. Figures 4.3 and 4.4 are two-dimensional and three-dimen- 
sional plots of the normal stress/maximum normal! stress for the five cases 
across the upper and lower interfaces, respectfully. For both interfaces, the 
maximum normal stress occurred in case 5 at the right hand side of the 
midlayer. Similar characteristics were displayed by all cases for both the 
upper and lower. As the thickness of the midlayer decreased, the ratio of the 
normal stress to the maximum normal stress decreased. Of interest was that 
the thinner material B became, the extent of the midspan that was at a 
minimum grew. Figures 4.5 and 4.6 show the bearing stresses acting along the 
upper and lower interfaces of material B. 

B. SHEAR STRESS DISTRIBUTION 

Figure 4.7 contains plots of the shear stresses across the upper and lower 
interfaces of material B for cases 5, 12, 13, 14, and 15. Cases 5 and 12 (with 
their greater thicknesses of material B) displayed similar attributes. The 
shear stresses, upper and lower, tended to axially compress the midlayer. 
Cases 13, 14, and 15 (with midlayer thicknesses of 0.25mm, 0.10mm, and 
0.05mm) showed peak shear stresses at the left and right ends of the 


interfaces, but along the majority of the length of span the shear stresses 
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Figure 4.3 2-D And 3-D Plots Of Normal Stress/Maximum Normal Stress At The 


Upper Interface For Cases 5, 12, 13, 14, And 15 
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Figure 4.4 2-D And 3-D Plots Of Normal Stress/Maximum Normal Stress At The 


Lower Interface For Cases 5, 12, 18, 14, And 15 
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fluctuated around zero. This fluctuation of positive and negative values of the 
shear stresses from point to point may have been due to a breakdown in the 
finite element program in dealing with the very small numerical values 
present. 

Comparisons of the ratio of the shear stress at a point to the maximum 
shear stress found in the group were conducted. Figures 4.8 and 4.9 are two- 
dimensional and three-dimensional plots of the shear stress/maximum shear 
stress for the five cases along the upper and lower interfaces, respectfully. For 
both the upper and lower, the maximum shear stress was located at the left 
end of material B. The two figures were similar, with the shear stress ratios 
being maximum at the extreme ends of the midlayer and then decreasing to 
a minimum towards midspan. Of interest, like the normal stress, was that as 
the thickness decreased, the range of the midspan that was at a minimum 
increased. Figures 4.5 and 4.6 display the shear stresses acting along the 
upper and lower interfaces for the midlayer. 

C. SUMMARY 

The effect of decreasing midlayer thickness on the tri-material configura- 
tion was displayed in a number of ways. The normal stresses remained 
bearing throughout as the midlayer became thinner. The maximum normal 
and shear stresses decreased with decreasing thickness. Minimal normal and 
shear stresses extended out further along the span and the stress distributions 


became more local in nature with the thinning of material B. A conclusion 
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Figure 4.8 2-D And 3-D Plots Of Shear Stress/Maximum Shear Stress At The 


Upper Interface For Cases 5, 12, 13, 14, And 15 
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Figure 4.9 2-D And 3-D Plots Of Shear Stress/Maximum Shear Stress At The 
Lower Interface For Cases 5, 12, 13, 14, And 15 


47 


that could be drawn would be that as the thickness of the midlayer continues 
to decrease, that the only stresses of any consequence along the span exist at 
the extreme left and right ends of the interfaces. The normal stresses were 
only one-third to one-half the shear stresses. For both the upper and lower 
interfaces the maximum normal stresses remained at the far left edge in all 
cases. The maximum shear stresses along the lower interface moved from the 
left edge to the right edge as the midlayer thickness was decreased. Finally, 
the maximum shear stresses along the upper interface remained at the left 


edge for all cases. 


48 


V. EFFECT OF MIDLAYER LENGTH 


In this chapter, the effect of the length of the midlayer (material B) on the 
tri-material configuration was investigated. Six different midlayer lengths 
were compared and studied. As shown in Figure 5.1, the following lengths for 
material B were used: 9mm (case 16), 7.5mm (case 17), 5mm (case 12), 2.5mm 
(case 18), 1mm (case 19), and 0.5mm (case 20). For all cases, the dimensions 
of materials A and C were kept constant (10mm length and 1mm thickness). 
Additionally, the thickness of the midlayer was held at 0.5mm and the tri- 
material configuration was given a temperature increase of 100°C. The 
material properties for the six cases were as follows: E,=E,=E,= 100 GPa, 
01,=100x10°/°C, a g= 300x10°/°C, and a = 10x10°/°C. 

A. NORMAL STRESS DISTRIBUTION 

Figure 5.2 contains plots of the normal stresses across the upper and 
lower interfaces of material B for the six cases. To allow comparison of the 
different cases the abscissa was constructed to be a non-dimensional location, 
with the distance along the interface being divided by the total midlayer length 
for that case. 

In all six cases, the normal stresses remained bearing across the entire 
span for both the upper and lower interfaces. Cases 12, 16, 17, and 18 (with 


midlayer lengths of 5mm, 9mm, 7.5mm, and 2.5mm, respectfully) displayed 
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Figure 5.1 Tri-material Configurations For Cases 12, 16, 17, 18, 19, And 20 
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similar characteristics. The maximum bearing stresses occurred at the left 
and right ends of the midlayer with the minimum being developed at midspan. 
For midlayer lengths of 9mm and 7.5mm, the stress along the midspan was 
negligible. As the length was further decreased to 5mm and 2.5mm, the 
midspan normal stresses increased. The stresses at the left and right ends 
simularly increased but to a much lesser extent as the length of the above 
mentioned four cases decreased. 

Cases 19 and 20, with the short midlayer lengths of Imm and 0.5mm, 
displayed distinctly different normal stress distributions than the four other 
cases. The left ends of the interfaces had relatively low bearing stresses. 
Across the remaining 90% of the span, the stress distribution was quite 
uniform. Figures 5.3 and 5.4 display the normal stresses acting along the 
upper and lower interfaces of the midlayer. 

B. SHEAR STRESS DISTRIBUTION 

Plots of the shear stresses along the upper and lower interfaces of the 
midlayer for cases 12, 16, 17, 18, 19, and 20 are displayed in Figure 5.5. Like 
the normal stress distribution plots, the abscissa was constructed to be a non- 
dimensional location, with the distance along the interface being divided by the 
total midlayer length for that case. 

All six cases showed similar attributes with both the upper and lower 
shear stresses tending to axially compress the midlayer. Additionally, the 


peak stresses occurred at the extreme left and right ends of the interfaces. For 
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Figure 5.5 Shear Stress At Interfaces As A Function Of Non-dimensional 
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the greatest length of material B, case 16 with L,= 9mm, the only non- 
negligible shear stresses appeared at the ends of the interfaces. As the 
midlayer length decreased, the extent of the span that experienced appreciable 
shear stress increased. This is clearly displayed in Figures 5.3 and 5.4. 
C. SUMMARY 

The effect of changing the midlayer length on the tri-material configura- 
tion was shown in numerous ways. As the material B length is shortened, the 
normal stresses tend to increase and become uniformly distributed across the 
span. The shear stress became distributed across a greater extent of the 
interface as the length of the midlayer lessened. The nature of the normal 
stresses was always bearing and were much less than the shear stresses 
generated. As the length of the midlayer was decreased to 1mm, shear and 
normal stresses increased. However, further reduction in length resulted in 
a slight decrease of both stresses. Except for the shortest length where it was 
at the right hand edge, the maximum shear stresses occurred at the left hand 
edge for the upper and lower interfaces. The maximum normal stresses on 
both interfaces stayed at the right edge. Finally, the maximum upper normal 
stresses were 10 to 30 percent greater than the maximum lower normal 


stresses. 
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VI. EFFECT OF NONUNIFORM TEMPERATURE DISTRIBUTION 


Chapters III, IV, and V considered the tri-material configuration as having 
a uniform temperature increase of 100°C (AT = 100°C). In this chapter the 
temperature increases of the three materials were allowed to vary, and the 
results were studied and compared. The physical dimensions of the tri- 
material configuration used in this chapter are shown in Figure 6.1. Materials 
A and C were both given a length of 10mm and a height of lmm. Material B 
was given a length of 2.5mm and a height of 0.1mm. The material properties 


supplied for all cases studied were as follows: E, = E, = E, = 100 GPa, a, = 


100x10°/°C, a = 300x10°/°C, and a, = 10x10°/°C. 
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Figure 6.1 Tyri-material Configuration For Cases 21, 22, 23, And 24 









Four different cases were examined. In case 21, all three materials were 
allowed a temperature rise of 100°C (AT, = AT, = AT, = 100°C). A tempera- 
ture increase of 200°C for materials A, B, and C (AT, = AT; = AT = 200°C) 


was studied in case 22. For case 23, material A was assigned a temperature 


oT 


rise of 200°C while materials B and C were permitted 100°C temperature 
increases (AT, = 200°C, AT; = AT, = 100°C). Finally, for case 24, material A 
was allowed to have a linearly decreasing temperature profile. The top layer 
of material A was set at a temperature increase of 200°C and the temperature 
rise was then lessened until at the interface of materials A and B the 
temperature increase was 100°C. This profile is shown in Figure 6.2. 


Materials B and C were allowed 100°C rises in temperature (AT, = linear, AT; 


= AT, = 100°C). 
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Figure 6.2 Material A Linear Temperature Distnibution For Case 24 
NORMAL STRESS DISTRIBUTION 


The normal stresses acting along the upper and lower interfaces of the 
midlayer for the four cases are shown in Figure 6.3. Similar characteristics 
were displayed for both the upper and lower normal stress distributions. The 
maximum stresses occurred at the far left and right ends of the interfaces and 
were always bearing. For the majority of the midspan the normal stresses 


were minimal. 
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Figure 6.3 Normal Stresses At The Interfaces For Cases 21, 22, 23, And 24 
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A comparison between case 21 and case 22 revealed an interesting point. 
Case 22, with a tri-material configuration temperature rise of 200°C, resulted 
in normal stresses which were twice the value of the normal stresses of case 
21 (AT, = AT, = AT, = 100°C). This showed the direct relationship between 
system temperature increase and the resulting increase in normal stresses. 

When only material A was given a larger temperature rise, as in case 23 
(AT, = 200°C, AT, = AT, = 100°C), the change in system behavior for normal 
stresses when compared toa uniform temperature field (case 21) was minimal. 
Along the upper interface the response was nearly identical, with the only 
difference being the marginally larger bearing stress at the extreme right end 
for case 23. Across the lower interface, case 23 normal stress values were 
continuously slightly less than those of case 21. 

When the temperature rise in material A was allowed to linearly decrease 
through its thickness, as in case 24, the resulting normal stress distributions 
were the lowest of the four cases studied. The maximum stresses developed 
were approximately one-half of the maximum for the uniform temperature field 
case. 

Figures 6.4 and 6.5 display the normal stress distributions acting on the 
upper and lower interfaces of the midlayer for the four cases. 

B. SHEAR STRESS DISTRIBUTION 
Figure 6.6 contains the plots of the shear stress distributions across the 


upper and lower interfaces for cases 21, 22, 23, and 24. The resulting curves 
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Figure 6.6 Shear Stresses At The Interfaces For Cases 21, 22, 23, And 24 







63 








2.5 


were similar for all four cases. The maximum shear stresses occurred at the 
extreme right and left ends of the interfaces. Along the majority of the length 
of span the shear stresses fluctuated around zero. This fluctuation may have 
been due to a breakdown in the finite element program in dealing with the 
very small numerical values present. 

As with the normal stresses, a comparison of case 21 (AT, = AT,;= AT, = 
100°C) with case 23 (AT, = AT, = AT, = 200°C) shows that the shear stresses 
at every point for the higher temperature increase were twice the value of 
those stresses for the lower temperature rise. This also displayed the direct 
relationship between system temperature increase and the resulting increase 
in shear stresses. 

Case 23 (AT, = 200°C, AT, = AT, = 100°C) and case 24 (AT, = linear, AT, 
= AT, = 100°C) produced shear stress distribution values very close to case 21. 
This suggests that for the particular tri-material configuration chosen for the 
chapter, increasing the temperature rise for one material does not appreciably 
change the system's shear stress response. 

Displays of the shear stress distributions acting on the upper and lower 
interfaces for the four cases are shown in Figures 6.4 and 6.5. 

C. SUMMARY 

A study of the effects of a nonuniform temperature field produced a 

number of interesting results. A doubling of the system temperature rise 


resulted in a doubling of the normal and shear stresses. Normal stresses 
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remained bearing in different uniform and nonuniform temperature fields. 
The normal and shear stress distributions maintained the same general shapes 
regardless of temperature distribution. Increasing the temperature rise in just 
one material resulted in minimal change in system response. Finally, a 
temperature rise which was allowed to linearly decrease through the thickness 


of one material resulted in the lowest normal stress distributions. 
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VII. SUMMARY AND CONCLUSIONS 


The results of this investigation into the behavior of a tri-material medium 
subjected to a temperature field leads to several conclusions. These 
conclusions are broken down into the following topics: effect of material 
properties, effect of thickness of midlayer, effect of length of midlayer, and 
effect of nonuniform temperature field. 

A. EFFECT OF MATERIAL PROPERTIES 

When both the temperature rise and the coefficients of thermal expansion 
were kept uniform throughout the system and only Young's Modulus was 
allowed to vary, no normal or shear stresses were developed. 

Changing the coefficient of thermal expansion produced several results. 
As the mismatch between the coefficients of thermal expansion ofthe materials 
became greater, the normal and shear stress magnitudes along the upper and 
lower interfaces became greater. The shear stresses developed along the 
interfaces were much greater than the normal stresses. Depending on the 
values of the coefficients of thermal expansion, the normal stresses were either 
bearing or peeling. Finally, along the upper and lower interfaces the normal 
and shear stresses were very local in nature. Only in the immediate area of 


the left and right corners of the interface were the stresses significant. 
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B. EFFECT OF MIDLAYER THICKNESS 

For the particular tri-material configuration studied the normal stresses 
remained bearing as the midlayer became thinner. The maximum normal and 
shear stresses decreased with decreasing material B thickness. As the 
midlayer was made thinner, the normal and shear stress distributions became 
more and more local in nature. It is believed that as the midlayer becomes 
extremely thin, the only stresses of any consequence would exist at the 
extreme left and right corners of the interfaces. The shear stresses generated 
were two to three times greater than the normal stresses. As the thickness of 
material B was decreased, the location of the maximum shear stress along the 
lower interface moved from the inner edge to the outer edge. However, along 
the upper interface the maximum shear stress remained fixed at the inner 
edge. Lastly, the maximum normal stress stayed at the left edge for both the 
upper and lower interfaces. 
C. EFFECT OF MIDLAYER LENGTH 

With the tri-material configuration used the resulting shear stresses were 
much greater than the normal stresses. The normal stresses developed were 
always bearing. As midlayer length was decreased to 1mm, the normal and 
shear stresses increased. Further shortening of the length of the midlayer 
resulted in slight decrease of the shear and normal stresses. On both 
interfaces the maximum shear stress occurred at the left hand edge. The only 


exception was the shortest length, where the maximum shear stress appeared 
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at the right hand edge. The maximum normal stress was experienced at the 
right hand edge for both the upper and lower interfaces. Additionally, the 
upper maximum normal stress was approximately 10 to 30 percent larger than 
the lower maximum normal stress. 
D. EFFECT OF NONUNIFORM TEMPERATURE DISTRIBUTION 

When the temperature rise for all three materials of the system was 
doubled, the result was a doubling of the normal and shear stresses. The 
normal stresses remained bearing and the shear and normal stress distribu- 
tions maintained the same general shapes regardless of the temperature 
distribution. Increasing the temperature rise in just the upper layer resulted 
in minimal change in Siem response. Finally, a linearly decreasing 
temperature rise gave the lowest magnitude normal! stress distributions. 
E. RECOMMENDATIONS FOR FURTHER RESEARCH 

A number of different areas for further research in this topic exist. 
Validate the tri-material response to nonuniform temperature fields by using 
a larger number of system configurations. Alter the finite element code to 
allow study of the system response to transient temperature fields. Refine the 
tri-material configuration to more closely match electronic component 


packaging. Finally, investigate the effect of thermal cycling on fatigue. 
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APPENDIX A 
FEM PROGRAM 


POUT SETS C TCC TSC CSCC CCSSSCCOSSCCSCCISCCOCSCSICSCS SOS SSCS SCSCSOCSCLCSS LASS eS eee ss 


me * 
* FEM PROGRAM FOR WELD PROBLEM : 
- Fa 


RAEARAKKKAAKRAAAEREKKKKAKKEKKKKaAKHEKRKEKKRAKKAKKEaEKAKKAKKEKKKKKKKKKKKKKKKkhkkkkkkkhkak 
* 


* COMPILE AND RUN THIS CODE WITH THE SINGLE EXEC COMMAND: WELD 
& 

& 

* NEL THE NUMBER OF ELEMENTS 

* NODES THE NUMBER OF SYSTEM NODES 

* NDOF THE NUMBER OF SYSTEM DEGREES OF FREEDOM 

* NBC THE NUMBER OF ESSENTIAL BOUNDARY CONDITIONS 
* IBC THE SUPPRESSED DEGREES OF FREEDOM 

* NFC THE NUMBER OF APPLIED FORCES 

* IFC THE SYSTEM NODES WITH APPLIED LOADS 

* ICORR THE CORRESPONDENCE BETWEEN LOCAL AND SYSTEM DOF 
*  NMAT NUMBER OF DIFFERENT MATERIALS 

*  IMAT A PARTICULAR MATERIAL 

* ITYPE = 0 FOR B.C.’S AT LEFT END ONLY, 

* = 1 FOR B.C.’S AT BOTH LEFT AND RIGHT ENDS 
* o£ YOUNGS MODULUS OF IMAT MATERIAL 

= G SHEAR MODULUS OF IMAT MATERIAL 

* XxX THE X COORDINATE OF A NODE 

* HITE THE HEIGHT OF AN ELEMENT 

*  ZLEN THE LENGTH OF AN ELEMENT 

* EK THE 6 BY 6 ELEMENT STIFNESS MATRIX 

* SYSMAT THE NDOF BY NDOF SYSTEM STIFFNESS MATRIX 

* FLOAD THE MAGNITUDE OF AN APPLIED LOAD 

* FORCE THE SYSTEM FORCE VECTOR 

* U THE SYSTEM DISPLACEMENT VECTOR 

* DELTMP THE TEMPERATURE CHANGE 

* ALPHA THERMAL COEFFICIENT OF EXPANSION 

* NNSCOL THE NUMBER OF NONSOLDERED COLUMNS 

* NSCOL THE NUMBER OF SOLDERED COLUMNS 

* N1ROW THE NUMBER OF ROWS IN THE CHIP 

* N3ROW THE NUMBER OF ROWS IN THE SUBSTRATE 

* PLSTRS THE PEELING STRESS AT THE NODE TO THE RIGHT OF THE 
* COLUMN NUMBER 

* N3STRT ELEMENT NUMBER BEFORE START OF LOWER SOLDER 
*& 


INCLUDE ‘WELDCOMN’ 
CHARACTER*20 NAME 
COMMON /WORKSP/RWKSP 
REAL RWKSP (137047) 
CALL IWKIN(137047) 


* 


READ *, NAME 


& 
* NOTE; INPUT THE NAME OF THE FILE TO BE OPENED AS FOLLOWS: 
* ‘/FN FT! 
* THE FILENAME MUST BE PRECEDED WITH A ’/ AND THE FILETYPE MUST BE 
* FOLLOWED BY A ’. THE ONLY SPACE IS BETWEEN FN AND FT. 
& 
OPEN(9, FILE= NAME ) 
* 
. INPUT THE FEM PARAMETERS 


READ(9,*) NEL, NODES, NMAT, NDOF, NBC, NFC, ITYPE 
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+ © 


+ 


‘ READ(9,*) (IBC(I), I = 1,NBC) 


DO 50 IB = 1,NBC 
50 CONTINUE 


IF (NFC .NE. 0) THEN 
READ(9,*) (IFC(I), FLOAD(I), I = 1,NFC) 
DO 70 IF = 1,NFC 
70 CONTINUE 
ELSE 
ENDIF 


DO 80 IM = 1,NMAT 

READ(9,81) E(IM), POISS(IM), ALPHA(IM), HITE(IM) 
80 CONTINUE 
81 FORMAT (4F10.7) 


READ(9,22) (X(IN), IN = 1,NODES) 
DO 999 ILL = 1, NEL 
READ(9,20) IL,IMAT(ILL), (ICORR(ILL,IN),IN = 1,6), INODE(ILL,1), 
& INODE(ILL, 2) 

999 CONTINUE 

20 FORMAT(10I4) 

22 FORMAT(7F10.5) 
READ(9,*) DELTMP 
READ(9,*) (TEMP(I),I = 1,NMAT+1) 
WRITE(14,*) ' ' 
WRITE(14,*) 'DELTMP = ’, DELTMP 


CALCULATE THE SHEAR MODULII 


WRITE(14,*%) ’ ’ 
WRITE(14,*) ’IMAT#, E, G MODULII, ALFA, POISSON AND HEIGHT FOLLOW’ 
WRITE(14,%) / ? 
DO 100 IM = 1, NMAT 
G(IM) = 5.*E(IM)/(12.*(1+POISS(IM) )) 
WRITE(14,105) IM, E(IM), G(IM), ALPHA(IM), POISS(IM), HITE(IM) 
100 CONTINUE 
105 FORMAT(3X, I3,2X,5E12.5) 


WRITE(14,*) ’ 
DO 350 I = 1,NDOF 
FORCE(I) = 0.0 
DO 340 J = 1,NDOF 
SYSMAT(I,J) = 0.0 
340 CONTINUE 
350 CONTINUE 
THE ELEMENT DO LOOP TO 500 FORMS THE SYSTEM MATRIX AND VECTOR 


DO 500 IEL = 1,NEL 
IM = IMAT(IEL) 


FORM THE ELEMENT STIFFNESS MATRICES 


IL = INODE(IEL, 1) 
IR = INODE(IEL, 2) 
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, 2LEN(IEL) = X(IR) - X(IL) 
H = HITE(IM) 


EE = E(IM) 

GG = G(IM) 

A = GG*ZLEN(IEL) / (4. *H) 
B = GG*H/ZLEN(IEL) 

C = EE*H/(6.*ZLEN(IEL) ) 
D = 0.5*GG 

AP2C =A + 2.%C 

ZMAPC = - A+C 

AM2C = A = 2.%C 

ZMAMC = - A- C 


EK(1,1) = AP2C 
EK(1,2) = ZMAPC 
EK(1,3) =D 
EK(1,4) = AM2C 
EK(1,5) = ZMAMC 
EK(1,6) = -D 
EK(2,1) = ZMAPC 
EK(2,2) = AP2C 
EK(2,3) = -D 
EK(2,4) = ZMAMC 
EK(2,5) = AM2C 
EK(2,6) = D 
EK(3,1) =D 
EK(3,2) = -D 
EK(3,3) = B 
EK(3,4) =D 
EK(3,5) = =D 
EK(3,6) = -B 
EK(4,1) = AM2C 
EK(4,2) = ZMAMC 
EK(4,3) =D 
EK(4,4) = AP2C 
EK(4,5) = ZMAPC 
EK(4,6) = =D 
EK(5,1) = ZMAMC 
EK(5,2) = AM2C 
EK(5,3) = -D 
EK(5,4) = ZMAPC 
EK(5,5) = AP2C 
EK(5,6) =D 
EK(6,1) = =D 
EK(6,2) = D 
EK(6,3) = -B 
EK(6,4) = -D 
EK(6,5) = D 
EK(6,6) = B 


* 


ASSEMBLE ELEMENT MATRICES INTO GLOBAL STIFFNES MATRIX 


DO 480 ID = 1,6 

IN = ICORR(IEL,ID) 

DO 480 JD = 1,6 

JN = ICORR(IEL,JD) 

SYSMAT(IN,JN) = SYSMAT(IN,JN) + EK(ID,JD) 
480 CONTINUE 
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(500 CONTINUE 
* 


* CONSTRUCT THE FORCE VECTOR DUE TO TEMPERATURE CHANGE 
® 
WRITE(14,*) ’ “ 
IF (DELTMP .NE. 0.0) THEN - 
DO 550 IEL = 1, NEL 
IM = IMAT(IEL) 
CHI2 = E(IM) *ALPHA(IM) *(2.*®TEMP (IM) +TEMP(IM+1) ) *HITE(IM) /6. 
CHI1 = E(IM) *ALPHA(IM) *(TEMP(IM) +2. *TEMP(IM+1) ) *HITE(IM) /6. 
I1= ICORR(IEL, 1) 
I2= ICORR(IEL, 2) 
I4= ICORR(IEL,4) 
I5* ICORR(IEL,5) 
FORCE(I1) = FORCE(I1) - CHI1 
FORCE(I2) = FORCE(I2) = CHI2 
FORCE(I4) = FORCE(I4) + CHI1 
FORCE(I5) = FORCE(IS) + CHI2 
550 CONTINUE 
ELSE 
ENDIF 


» 


MODIFY STIFFNESS MATRIX AND FORCE VECTOR FOR BOUNDARY CONDITIONS 


DO 610 I = 1,NBC 
II = IBC(I) 
DO 600 J = 1,NDOF | 
SYSMAT(II,J) = 0.0 
600 CONTINUE 
SYSMAT(II,II) = 1. 
FORCE(II) = 0.0 
610 CONTINUE 


IF(ITYPE .EQ. 1) THEN 
READ(9,*) ILASTTOP 
WRITE(14,*) ’ “ 
WRITE(14,*) ‘B.C. S AT BOTH ENDS; ILASTTOP = ’, ILASTTOP 
DO 620 IN = 1,NDOF 
SYSMAT(ILASTTOP,IN) = 0.0 
620 CONTINUE 
SYSMAT(ILASTTOP,ILASTTOP) = 1. 
SYSMAT(ILASTTOP,NDOF) = -1. 
FORCE(ILASTTOP) = 0.0 
ELSE 
ENDIF 


IF (NFC .NE. 0) THEN 
DO 630 IF = 1,NFC 
II = IFC(IF) 
FORCE(II) = FORCE(II) + FLOAD(IF) 
630 CONTINUE 
ELSE 
ENDIF 


» 


SOLVE THE SYSTEM OF EQUATIONS FOR THE DISPLACEMENTS 


CALL DLSARG(NDOF, SYSMAT, NPARAM, FORCE, 1, UV) 


* 


PRINT THE FORCE AND DISPLACEMENT VECTORS 
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WRITE(14,*) ’ “ 
, WRITE(14,*) ‘ FORCE AND DISPLACEMENT VECTORS FOLLOW’ 
WRITE(14,*) ’ ' 
DO 680 IDOF = 1,NDOF 
WRITE(14,641) IDOF, FORCE(IDOF), U(IDOF) 
680 CONTINUE 
641 FORMAT(2X, 14, 4X, E12.5, 10X, E12.5) 
x 
* CALCULATE THE STRESSES 
« 
DO 700 IEL = 1, NEL 
I1 = ICORR(IEL, 1) 
ICORR(IEL, 2) 
ICORR(IEL, 3) 
ICORR(IEL, 4) 
ICORR(IEL,5) 
ICORR (IEL, 6) 
U1BOT = U(I1) 
U1TOP = U(I2) 
V1 = U(I3) 
U2BOT = U(I4) 
U2TOP = U(IS5) 
v2 = U(I6) 
IM = IMAT(IEL) 
EE = E(IM) 
GG = G(IM) 
SIGMABOT(IEL)= EE*((U2BOT - U1BOT)/ZLEN(IEL) - ALPHA(IM) * 
& TEMP(IM+1) ) 
SIGMATOP(IEL)= EE*((U2TOP - U1TOP)/ZLEN(IEL) - ALPHA(IM) *TEMP (IM) ) 
TAUELEM(IEL) = GG*( ((U1TOP -U1BOT) + (U2TOP - U2BOT) ) / 
& (2.*HITE(IM)) + (V2 - V1)/ZLEN(IEL) ) 
700 CONTINUE 


HHH 
mwWn 
hann 


i 
ov 
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WRITE(14,*) ’ / 

WRITE(14,*) “ STRESSES FOLLOW’ 

WRITE(14,*) ’ 

WRITE(14,*) / SGMABOT- - SGMATOP TAUELEM ’ 


PSI = E(1) *ALPHA(1) *DELTMP 

DO 800 IEL = 1,NEL 

IM = IMAT(IEL) 

SIGBTNDL = SIGMABOT(IEL) /PSI 

SIGTPNDL = SIGMATOP(IEL) /PSI 

TAUNDL = TAUELEM(IEL) /PSI 

WRITE(14,801) IEL, SIGMABOT(IEL), SIGMATOP(IEL), TAUELEM(IEL) 


800 CONTINUE 
801 FORMAT ((2X,14,6(4X,G12.5))) 


WRITE(14,*) ’ / 

READ(9,*) NNSCOL, NSCOL, N1ROW, N3ROW 
NCOL = NNSCOL + NSCOL 

DO 1300 IEL = 6, 70 

IM = IMAT(IEL) 

PLFCE(IEL) = TAUELEM(IEL) *HITE(IM) 


1300 CONTINUE 
II = 6 
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3 
* 


* 
* 


c 
* 


1400 


1600 


1700 


1800 
2 


DO 1500 ICOL = 1,8 

IELF = 6 + (ICOL = 1) 

IELL = IELF + 64 

PLF = 0.0 

DO 1400 IEL = IELF,IELL, 14 

PLF = PLF + PLFCE(IEL+1) = PLFCE(IEL) 
CONTINUE 

UPLSTR(ICOL) = 2.*PLF/(X(II+2) = X(II)) 
WRITE(14,31) ICOL, UPLSTR(ICOL) 

Il =e + 

CONTINUE 

FORMAT (1X,I2,’ UPPER PEEL STRESS = ’,G12.5) 


WRITE(14,%) ’% % 

DO 1600 IEL ® 121, 197 

IM = IMAT(IEL) 

PLFCE(IEL) = TAUELEM(IEL) *HITE(IM) 

CONTINUE 

II = 6 

DO 1800 ICOL = 1,8 

IELF = 121 + (ICOL = 1) 

IELL = IELF + 76 

PLF = 0.0 

DO 1700 IEL = IELF,IELL,17 

PLF = PLF + PLFCE(IEL+1) - PLFCE(IEL) 
CONTINUE 

BPLSTR(ICOL) #(2.*®PLF/(X(II+2) = X(II)))*(-1.0) 
WRITE(14,32) ICOL, BPLSTR(ICOL) 

Il "=a + 1 

CONTINUE 

FORMAT (1X,I2,’ LOWER PEEL STRESS = ’,G12.5) 


COMPARE MOMENTS DUE TO BENDING AND PEEL STRESSES 


1900 


2000 


AVEFCE(1) = 0.5*(SIGMATOP(1) + SIGMABOT(1)) *HITE(1) 
ZMOMARM(1) = .0562 - .014_ 

AVEFCE(2) = 0.5*(SIGMATOP(15) + SIGMABOT(15) ) *HITE(2) 
ZMOMARM(2) = .0562 - .028 - .007 

AVEFCE(3) = 0.5*(SIGMATOP(29) + SIGMABOT(29) ) *HITE(3) 
ZMOMARM(3) = .0562 - .028 - .014 = .0035 

AVEFCE(4) = 0.5*(SIGMATOP(43) + SIGMABOT(43) ) *HITE(4) 
ZMOMARM(4) = .0035 + .00175 

AVEFCE(5) = 0.5*(SIGMATOP(57) + SIGMABOT(57) ) *HITE(4) 
ZMOMARM(5) = .00175 

ZMOMBDG = 0.0 

DO 1900 I = 1,5 

ZMOMBDG = ZMOMBDG + ZMOMARM(I) *AVEFCE(I) 

CONTINUE 


ZMOMPEEL = 0.0 
DO 2000 I = 1,8 


II =7+HtI 
JJ = il - 2 
KK = II - 1 


W= 0.5*(X(II) - X(JJ)) 

Y = 0.25*(X(JJ) + X(II) +2.*X(KK) ) 
ZMOMPEEL = ZMOMPEEL + W*Y*UPLSTR(TI) 
CONTINUE 

WRITE(14,*%) ‘MOMPEEL = ’, ZMOMPEEL 
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(ZMOMBDG = 0.0 
IM = 0 

Y2 = 0.056 

DO 3000 I = 1,57,14 

IM = IMAT(I) 

Yl = Y2 - HITE(IM) 

AOV3 = (SIGMATOP(I)-SIGMABOT(I))/(3.*HITE(IM) ) 

BOV2 = 0.5*(SIGMATOP(I)-(Y¥2/HITE(IM) ) * (SIGMATOP(I) -SIGMABOT(I)) ) 
ZMOMBDG = ZMOMBDG + AOV3*(Y2**3-¥1**3) +BOV2*(Y2**2-Y1"*2) 

Y2 = Yl 

3000 CONTINUE 

SIGMAT = 1.052 

SIGMAB = 2.3012 

Y2 = .001 

Yl = .0005 

B = SIGMAT - (Y2/.0005) *(SIGMAT-SIGMAB) 

A = (SIGMAT-SIGMAB) /(.0005) 

STOP 

END 
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APPENDIX B 
PROGRAM EXEC 


FORTVS2 WELD (AD(DBL4) ICA) 
FILEDEF 14 DISK WELD OUTPUT (RECFM FB LRECL 102 BLOCK 10200) 
LOAD WELD (START 
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APPENDIX C 
INPUT FILE 
200 18 15 261 14 0 0 


123 45 6 7 146 164 165 166 167 168 169 
100.E+09 .3 100.E-+06 ~500 












100.E+09 -3 100.E-06 ~250 
100.E+09 -3 100.E-06 25 
100.E+09 -3 100.E-06 -0625 
100.E+09 -3 100.E-06 -0625 
100.E+09 -3 300.E-06 -010 
100.E+09 -3 300.E-06 ~020 
100.E+09 2 -3 300.E-06 3 .040 4 


100.E+09 -3 300.E-06 -020 


100.E+09 .3 300.E-06 010 
100.E+09 .3. 10.E-06 0625 
100.E+09 .3 »=« 10. E-06 0625 
100.E+09 .3. 10.E-06 £125 
100. E+09 .3 10sE-06 .250 
00.E+09 -3 10.E-06, 500, 
0.0 1.33333 2.66667 4.0 4.66667 5.0 Sos 
Ba50 6.0 7.0 8.0 9.0 9.50 9.750 
10. 10.000001 10.000002 10.000003 
7) Se 1m 6 95 $= @ 2 
> @ We Ge yomeiee 15 16mm > ° 3 
A 91 tee 1S. temeee os 23. Bo 4 
4 1 24 22 23 31 29 30 4 «5 
5 1 31 29 30 38 36 37 5 6 
6 1 38 36 37 45 43 44 #6 7 
7 1 45 43 44 S2 50 51 7 8 
8 1 52 50 51 59 57 58 8 9 1. Young’s Modulus 
6 4 6S S756 6a Gufe65 9 10 
10 1 66) 6265 F307? we Go 11 ae 
11 1 73 #71 «#72 «#80 78 #79 #11 © «12 2. Poisson's Ratio 
12 1 80 78 79 87 85 86 12 13 
13. 1 Gtme6566)8Oh S208os 43 Te BeGceticent of 
14 1 94 92 93 101 99100 14 15 . 
15 2 eo Mi io es 1 == 2 Thermal Expansion 
16 2 11810 WO Pielody) 16 2 es 
17 2 #18 17 16 25 2% 23 3 4 4, Row Height 
18 2 25 4% S80 320uee 30 &@ Ss 
19 2 32 31 30 39 38 37 #5 6 
20 2 39 38 37 46 45 44 #6 7 5. Column Width 
21 #2 46 45 @4 53 82 51 7 8 
22 2 53 52 51 60 59 58 8 9 
23 2 60 59 58 67 66 65 9 10 
24 2 67 66 65 74 73 72 #10 11 
25 2 74 73 72 81 80 79 11 «12 
26 2 81 80 79 88 87 86 12 13 
27 2 88 87 86 95 94 93 13 14 
28 2 95 94 93 102 101 100 14 15 
a ee SW 8 ee 
$043 «+@ wi 9 19 16 16 2 3 
S143 «fees 16°26 25 33 3 4 
32 3 26 25 23 33 32 30 Sa 5 
33 3 33 32 30840 39 37 #5 6 
34 3 40 39 37 #47 46 «44«~«6€6~«=27 
35 3 WT veé wee WES Se “Sa 7 8 
36 3 SS0ss wa Gis GQ, 5S 4868 3 
37. 3 61860 56 66 6% 65 9 10 
98 3 68067 65 75 74) 72 10 12 
39 3 #75 #74 #72 #82 81 #79 #11 «12 
40 3 82 81 79 89 88 86 12 13 


ua 


OWU WD OWOODWOODODOOAINIIVNVI IY IWVIAHAAHAAAHAAAMNYNMNYNUNUNUNHNNNUUNUN MNS SG HS oS HHH FS HLL LLLwUWwY 


78 


79 


KH OWwWWO WAN AWN 


he 


161 13 226 225 79 231 230 eGee ee 
162‘ 13 231 230 86 236 235 93 13 14 
163 13 236 235 93 241 240 100 14 15 
164 13 241 240 100 247 246 244 15 16 


165 13 247 246 244 253 252 250 16 17 
166 13 253 252 250 259 258 256 17 18 
167 14 168 167 164 174 173 170 ut 2 
168 14 174 173 170 180 179 176 2 3 
169 14 180 179 176 186 185 182 8 4 
170 14 186 185 182 192 191 188 4 5 
171 #14 192 191 188 197 196 37 > 6 
172 14 197 196 37 202 201 44 6 7 
1739 14 202 201 44 207 206 51 7 8 
174 14 207 206 51 212 211 #458 8 9 
175 14 212 211 5S6 217 216 65 9 10 
176 14 217 216 65 222 221 72 10 «121 


177 14 222 221 72 227 226 FO ii )i2 
178 14 227 226 79 232 231 86 12 13. 
179 14 232 231 86 237 236 93 13 14 


180 14 237 236 93 242 241 100 14 15 
181 14 242 241 100 248 247 244 15 16 
182 14 248 247 244 254 253 250 16 17 
183 14 254 253 250 260 259 256 17 18 
184 15 169 168 164 175 174 170 1 2 
185 15 175 174 170 181 180 176 3 
186 15 181 180 176 187 186 182 3 4 
187 15 187 186 182 193 192 188 4 5 
188 15 193 192 188 198 197 37 5 6 
189 15 198 197 37 203 202 44 6 7 
190 15 203 202 44 208 207 51 ri 8 
191 15 208 207 Si 273 212556 8 9 
192 15 213 212 .@6 218 277 cs 9 10 
193. 15 218 217 65 223 222 72 10 11 
194 15 223 222 72 228 227 579i 
195 15 228 227 79 233 232) S62 lo 


196 15 233 232 86 238 237 93 13 14 

197 15 238 237 93 243 242 100 14 15 

198 15 243 242 100 249 248 244 15 16 

199 15 249 248 244 255 254 250 16 17 

200 15 255 254 250 261 260 256 17 18 

200. | 
100. 100. 100. 100. 100. 100. 100. 100. 100. 100. 100. 100. 100. —— 
100. 100. 
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6. Layer Temperature Rise 
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DELTMP #@ 


IMAT#, E 


APPENDIX D 
OUTPUT FILE 


200.000000000000000 


, G MODULII, ALFA, POISSON AND HEIGHT FOLLOW 


0.10000E+12 
0.10000E+12 
0.10000E+12 
0.10000E+12 
0.10000E+12 
0.10000E+12 
0.10000E+12 
0.10000E+12 
0.10000E+12 


0.10000E+12 


0.10000E+12 


0.10000E+12 


0.10000E+12 


0.10000E+12 


0.10000E+12 


0.32051E+11 
0.32051E+11 
0.32051E+11 
0.32051E+11 
0.32051E+11 
0.32051E+11 
0.32051E+11 
0.32051E+11 
0.32051E+11 
0.32051E+11 
0.32051E+11 
0.32051E+11 
0.32051E+11 
0.32051E+11 


0.32051E+11 


FORCE AND DISPLACEMENT VECTORS 


0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 


0.00000E+00 


81 


0.10000E-03 
0.10000E-03 
0.10000E-03 
0.10000E~-03 
0.10000E~-03 
0.30000E-03 
0.30000E-03 
0.30000E-03 
0.30000E-03 
0.30000E-03 
0.10000E-04 
0.10000E-04 
0.10000E-04 
0.10000E-04 


0.10000E-04 


FOLLOW 


0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 


0.00000E+00 


0.30000E+00 


0.30000E+00 


0.30000E+00 


0.30000E+00 


0.30000E+00 


0.30000E+00 


0.30000E+00 


0.30000E+00 


0.30000E+00 


0.30000E+00 


0.30000E+00 


0.30000E+00 


0.30000E+00 


0.30000E+00 


0.30000E+00 


0.50000E+00 


0.25000E+00 


0.12500E+00 © 


0.62500E-01 
0.62500E-01 
0.10000E-01 
0.20000E-01 
0.40000E-01 
0.20000E-01 
0.10000E-01 
0.62500E-01 
0.62500E-01 
0.12500E+00 
0.25000E+00 


0.50000E+00 


0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0. 00000E+00 
0.00000E+00 
0. 00000E+00 
0.Q00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0. 00000E+00 
0.00000E+00 
0.00000E+00 
0. 00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0. 00000E+00 


0.00000E+00 
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0.00000E+00 
0.00000E+00 


0.16022E-01 


-0.51672E-02 


0.12042E-01 
0.10033E-01 
0.90521E-02 
0.85658E-02 
0.80331E-02 


0.32451E-01 


-0.21458E-01 


0.24079E-01 
0.19888E-01 
0.17771E-01 
0.16716E-01 
0.15765E-01 


0.49296E-01 


-0.50058E-01 


0.36110E-01 
0.29507E-01 
0.26255E-01 
0.24626E-01 
0.22845E-01 


0.57798E-01 


-0.69213E-01 


0.42198E-01 
0.34279E-01 
0.30238E-01 
0.28255E-01 


0.26546E-01 


36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 


65 


0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0-00000E+00 
-0.15000E+08 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.Q00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.Q00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0. 00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 


0.00000E+00 


83 


0.61934E-01 
-0.80181E-01 
0.45342E-01 
0.36811E-01 
0.32342E-01 
0.29943E-01 
0.27144E-01 
0.65005E-01 
-0.88399E-01 
0.47763E-01 
0.38867E-01 
0.34269E-01 
0.31904E-01 
0.29455E-01 
0.68161E-01 
-0.97019E-01 
0.50215E-01 
0.41041E-01 
0.36432E-01 
0.34170E-01 
0.32225E-01 
0.74620E-01 
-0.11546E+00 
0.55219E-01 
0.45459E-01 
0.40623E-01 
0.38212E-01 
0.35508E-01 
0.87846E-01 


-0.15715E+00 
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67 


68 


69 


70 


va 


V2 


gs 


74 


75 


76 


77 


78 


79 


80 


81 


82 


83 


84 


85 


86 


87 


88 


89 


90 


91 


92 


93 


94 


95 


0.00000E+00 
0. 00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.Q00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.Q00000E+00 
0.Q00000E+00 
0.00000E+00 
0.Q00000E+00 
0.00000E+00 
0.Q00000E+00 
0.Q00000E+00 
0. 00000E+00 
0. 00000E+00 
0.00000E+00 
0.Q00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 


0.00000E+00 


0.65382E-01 
0.54158E-01 
0.48522E-01 
0.45704E-01 
0.43179E-01 


0.10123E+00 


-0.20527E+00 


0.75574E-01 
0.62744E-01 
0.56355E-01 
0.53160E-01 
0.49681E-01 


0.11433E+00 


-0.25976E+00 


0.85778E-01 
0.71538E-01 
0.64345E-01 
0.60737E-01 
0.57389E-01 


0.12051E+00 


-0.28934E+00 


0.90866E-01 
0.76222E-01 
0.68790E-01 
0.64988E-01 
0.60910E-01 


0.12339E+00 


-0.30467E+00 


0.93387E-01 


0.78682E-01 
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96 
97 
98 
99 
100 
101 
102 
103 
104 
105 
106 
107 
108 
109 
110 
111 
112 
113 
114 
175 
116 
117 
118 
119 
120 
a ap | 
122 
123 
124 


125 


0.00000E+00 
0.00000E+00 
0.00000E+00 
0.25000E+09 
0.00000E+00 
0.37500E+09 
0.18750E+09 
0.93750E+08 
0.62500E+08 
0.46250E+08 
-0.45000E+08 
0.00000E+00 
0.00000E+00 
0. OOOO0E+00 
0.00000E+00 
0. OOO00E+00 
0.00000E+00 
0.00000E+00 
0.Q00000E+00 
0.45000E+08 
-0.90000E+08 
0. 00000E+00 
0. 00000E+00 
0.00000E+00 
0. 00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 


0.90000E+08 
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0.71308E-01 
0.67581E-01 
0.63741E-01 
0.12604E+00 
-0.32032E+00 
0.95892E-01 
0.81178E-01 
0.73861E-01 
0.70313E-01 
0.67120E-01 
0.25831E-01 
0.29799E-01 
0.31288E-01 
0.35675E-01 
0.42187E-01 
0.49705E-01 
0.56277E-01 
0.60826E-01 
0.62514E-01 
0.67226E-01 
0.24860E-01 
0.29050E-01 
0.30760E-01 
0.34696E-01 
0.41504E-01 
0.48460E-01 
0.55339E-01 
0.59371E-01 
0.61365E-01 


0.65955E-01 


126 
127 
128 
129 
130 
32 
132 
133 
134 
135 
136 
137 
138 
139 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
151 
152 
153 
154 


155 


-0.90000E+08 


0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 


0.90000E+08 


-0.45000E+08 


0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.45000E+08 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 


0.00000E+00 


-0.15000E+08 


0.00000E+00 


0.00000E+00 


0.00000E+00 


0.00000E+00 
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0.23563E-01 
0.27710E-01 
0.29341E-01 
0.33167E-01 
0.39696E-01 
0.46422E-01 
0.53003E-01 
0.56897E-01 
0.58707E-01 
0.63070E-01 
0.23197E-01 
0.27159E-01 
0.28411E-01 
0.32652E-01 
0.38538E-01 
0.45660E-01 
0.51577E-01 
0.55900E-01 
0.57197E-01 
0.61436E-01 
0.00000E+00 
0.71778E-02 
0.14053E-01 
0.20279E-01 
0.23548E-01 
0.23941E-01 
0.26069E-01 
0.28720E-01 
0.31633E-01 


0.38716E-01 


156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 
171 
172 
a7 3 
174 
175 
176 
177 
178 
179 
180 
181 
182 
183 
184 


185 


0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.15000E+08 
0.Q00000E+00 
0.00000E+00 
0.31250E+07 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.Q00000E+00 
0.00000E+00 
0.Q00000E+00 
0.00000E+00 
0.000Q00E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 


0.00000E+00 
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0.44524E-01 
0.51618E-01 
0.54641E-01 
0.57231E-01 
0.59728E-01 
0.59728E-01 
0.59728E-01 
0.59728E-01 
0.00000E+00 
0.00000E+00 
0.Q00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
-0.62365E-02 
0.66407E-02 
0.60583E-02 
0.49010E-02 
0.26339E-02 
=O. 19392E-02 
-0.24155E-01 
0.12867E-01 
O211783E=-01 
0.96238E-02 
0.52627E-02 
-0.34704E-02 
-0.52573E-01 
0.18851E-01 
0.17274E-01 


0.14111E-01 


186 
187 
188 
189 
190 
191 
192 
193 
194 
igo 
196 
197 
198 
199 
200 
201 
202 
203 
204 
205 
206 
207 
208 
209 
210 
211 
212 
213 
214 


2r5 


0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 
0.00000E+00 


0.00000E+00 


0.78855E-02 


-0.45865E-02 


-0.70479E-01 
0.21520E-01 
0.19760E-01 
0.16317E-01 
0.92685E-02 

-0.50653E-02 
0.22721E-01 
0.21115E-01 
0.17566E-01 
0.10061E-01 

-0.54197E-02 
0.24327E-01 
0.22481E-01 
0.18660E-01 
0.10718E-01 

-0.57161E-02 
0.26229E-01 
0.24082E-01 
0.19872E-01 
0.11406E-01 

-0.59275E-02 
0.29553E-01 
0.27155E-01 
0.22365E-01 
0.12883E-01 

-0.62042E-02 
0.35599E-01 


0.32803E-01 


88 


216 
217 
218 
219 
220 
221 
222 
223 
224 
225 
226 
227 
228 
229 
230 
231 
232 
233 
234 
235 
236 
237 
238 
239 
240 
241 
242 
243 
244 


245 


0. 00000E+00 
0. 00000E+00 
0. 00000E+00 
0.00000E+00 
0.00000E+00 
0. 00000E+00 
0.00000E+00 
0.00000E+00 
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